A Genome-Wide Arrayed CRISPR Screen Reveals PLSCR1 as an Intrinsic Barrier to SARS-CoV-2 Entry (preprint)

Interferons (IFNs) play a crucial role in the regulation and evolution of host-virus interactions. Here, we conducted a genome-wide arrayed CRISPR knockout screen in the presence and absence of IFN to identify human genes that influence SARS-CoV-2 infection. We then performed an integrated analysis of genes interacting with SARS-CoV-2, drawing from a selection of 67 large-scale studies, including our own. We identified 28 genes of high relevance in both human genetic studies of COVID-19 patients and functional genetic screens in cell culture, with many related to the IFN pathway. Among these was the IFN-stimulated gene PLSCR1. PLSCR1 did not require IFN induction to restrict SARS-CoV-2 and did not contribute to IFN signaling. Instead, PLSCR1 specifically restricted spike-mediated SARS-CoV-2 entry. The PLSCR1-mediated restriction was alleviated by TMPRSS2 over-expression, suggesting that PLSCR1 primarily restricts the endocytic entry route. In addition, recent SARS-CoV-2 variants have adapted to circumvent the PLSCR1 barrier via currently undetermined mechanisms. Our study contributes to understanding the association between PLSCR1 variants and severe COVID-19 cases reported in a recent GWAS.

The risk of COVID-19 death is much greater and age-dependent with type I IFN autoantibodies (preprint)

SARS-CoV-2 infection fatality rate (IFR) doubles with every five years of age from childhood onward. Circulating autoantibodies neutralizing IFN-α, IFN-ω, and/or IFN-β are found in ~20% of deceased patients across age groups. In the general population, they are found in ~1% of individuals aged 20-70 years and in >4% of those >70 years old. With a sample of 1,261 deceased patients and 34,159 uninfected individuals, we estimated both IFR and relative risk of death (RRD) across age groups for individuals carrying autoantibodies neutralizing type I IFNs, relative to non-carriers. For autoantibodies neutralizing IFN-α2 or IFN-ω, the RRD was 17.0[95% CI:11.7-24.7] for individuals under 70 years old and 5.8[4.5-7.4] for individuals aged 70 and over, whereas, for autoantibodies neutralizing both molecules, the RRD was 188.3[44.8-774.4] and 7.2[5.0-10.3], respectively. IFRs increased with age, from 0.17%[0.12-0.31] for individuals <40 years old to 26.7%[20.3-35.2] for those ≥80 years old for autoantibodies neutralizing IFN-α2 or IFN-ω, and from 0.84%[0.31-8.28] to 40.5%[27.82-61.20] for the same two age groups, for autoantibodies neutralizing both molecules. Autoantibodies against type I IFNs increase IFRs, and are associated with high RRDs, particularly those neutralizing both IFN-α2 and -ω. Remarkably, IFR increases with age, whereas RRD decreases with age. Autoimmunity to type I IFNs appears to be second only to age among common predictors of COVID-19 death.

Resolving the Functional Significance of BRCA1 RING Domain Missense Substitutions (preprint)

Part 1: Development and calibration of suitably accurate functional assays for BRCA1 RING domain and BRCT domain missense substitutions could dramatically accelerate clinical classification of rare missense substitutions observed in that gene. Leveraging data from 68,000 full sequence tests of BRCA1 and BRCA2, plus data from the limited number of already classified BRCA1 RING domain missense substitutions, we used logistic regression and related techniques to evaluate three BRCA1 RING domain assays. These were recently described high throughput yeast 2-hybrid and E3 ubiquitin ligase assays, plus a newly developed mammalian 2- hybrid assay. While there were concerns about the accuracy of the yeast 2-hybrid assay and the indirect nature of the ubiquitin ligase assay, the mammalian 2-hybrid assay had excellent correlation with existing missense substitution classifications. After calibration, this assay contributed to classification of one newly reported BRCA1 missense substitution. In principal, the mammalian 2-hybrid assay could be converted to a high-throughput format that would likely retain suitable accuracy. Part 2: How does one achieve clinically applicable classification of the vast majority of all possible sequence variants in disease susceptibility genes? BRCA1 is a high-risk susceptibility gene for breast and ovarian cancer. Pathogenic protein truncating variants are scattered across the open reading frame, but all known missense substitutions that are pathogenic because of missense dysfunction are located in either the amino-terminal RING domain or the carboxy-terminal BRCT domain. Heterodimerization of the BRCA1 and BARD1 RING domains is a molecularly defined obligate activity. Hence, we tested every BRCA1 RING domain missense substitution that can be created by a single nucleotide change for heterodimerization with BARD1 in a Mammalian 2-hybrid (M2H) assay. Downstream of the M2H laboratory assay, we addressed three additional challenges: assay calibration, validation thereof, and integration of the calibrated results with other available data such as computational evidence and patient/population observational data to achieve clinically applicable classification. Overall, we found that about 20% of BRCA1 RING domain missense substitutions are pathogenic. Using a Bayesian point system for data integration and variant classification, we achieved clinical classification of about 89% of observed missense substitutions. Moreover, among missense substitutions not present in the human observational data used here, we find an additional 47 with concordant computational and functional assay evidence in favor of pathogenicity; these are particularly likely to be classified as Likely Pathogenic once human observational data become available.

DRUL for School: Opening Pre-K with safe, simple, sensitive saliva testing for SARS-CoV-2 (preprint)

To address the need for simple, safe, sensitive, and scalable SARS-CoV-2 tests, we validated and implemented a PCR test that uses a saliva collection kit use at home. Individuals self-collected 300 ul saliva in vials containing Darnell Rockefeller University Laboratory (DRUL) buffer and extracted RNA was assayed by RT-PCR (the DRUL saliva assay). The limit of detection was confirmed to be 1 viral copy/ul in 20 of 20 replicate extractions. Viral RNA was stable in DRUL buffer at room temperature up to seven days after sample collection, and safety studies demonstrated that DRUL buffer immediately inactivated virus at concentrations up to 2.75x106 PFU/ml. Results from SARS-CoV-2 positive nasopharyngeal (NP) swab samples collected in viral transport media and assayed with a standard FDA Emergency Use Authorization (EUA) test were highly correlated with samples placed in DRUL buffer. Direct comparison of results from 162 individuals tested by FDA EUA oropharyngeal (OP) or NP swabs with co-collected saliva samples identified four otherwise unidentified positive cases in DRUL buffer. Over six months, we collected 3,724 samples from individuals ranging from 3 months to 92 years of age. This included collecting weekly samples over 10 weeks from teachers, children, and parents from a pre-school program, which allowed its safe reopening while at-risk pods were quarantined. In sum, we validated a simple, sensitive, stable, and safe PCR-based test using a self-collected saliva sample as a valuable tool for clinical diagnosis and screening at workplaces and schools.

Genome-scale identification of SARS-CoV-2 and pan-coronavirus host factor networks (preprint)

The COVID-19 pandemic has claimed the lives of more than one million people worldwide. The causative agent, SARS-CoV-2, is a member of the Coronaviridae family, which are viruses that cause respiratory infections of varying severity. The cellular host factors and pathways co-opted by SARS-CoV-2 and other coronaviruses in the execution of their life cycles remain ill-defined. To develop an extensive compendium of host factors required for infection by SARS-CoV-2 and three seasonal coronaviruses (HCoV-OC43, HCoV-NL63, and HCoV-229E), we performed parallel genome-scale CRISPR knockout screens. These screens uncovered multiple host factors and pathways with pan-coronavirus and virus-specific functional roles, including major dependency on glycosaminoglycan biosynthesis, SREBP signaling, and glycosylphosphatidylinositol biosynthesis, as well as an unexpected requirement for several poorly characterized proteins. We identified an absolute requirement for the VTT-domain containing protein TMEM41B for infection by SARS-CoV-2 and all other coronaviruses. This human Coronaviridae host factor compendium represents a rich resource to develop new therapeutic strategies for acute COVID-19 and potential future coronavirus spillover events. HIGHLIGHTSGenome-wide CRISPR screens for SARS-CoV-2, HCoV-OC43, HCoV-NL63, and HCoV-229E coronavirus host factors. Parallel genome-wide CRISPR screening uncovered host factors and pathways with pan-coronavirus and virus-specific functional roles. Coronaviruses co-opt multiple biological pathways, including glycosaminoglycan biosynthesis, SREBP signaling, and glycosylphosphatidylinositol biosynthesis and anchoring, among others. TMEM41B - a poorly understood factor with roles in autophagy and lipid mobilization - is a critical pan-coronavirus host factor.

Functional interrogation of a SARS-CoV-2 host protein interactome identifies unique and shared coronavirus host factors (preprint)

The ongoing SARS-CoV-2 pandemic has devastated the global economy and claimed nearly one million lives, presenting an urgent global health crisis. To identify host factors required for infection by SARS-CoV-2 and seasonal coronaviruses, we designed a focused high-coverage CRISPR-Cas9 library targeting 332 members of a recently published SARS-CoV-2 protein interactome. We leveraged the compact nature of this library to systematically screen four related coronaviruses (HCoV-229E, HCoV-NL63, HCoV-OC43 and SARS-CoV-2) at two physiologically relevant temperatures (33 {degrees}C and 37 {degrees}C), allowing us to probe this interactome at a much higher resolution relative to genome scale studies. This approach yielded several new insights, including unexpected virus and temperature specific differences in Rab GTPase requirements and GPI anchor biosynthesis, as well as identification of multiple pan-coronavirus factors involved in cholesterol homeostasis. This coronavirus essentiality catalog could inform ongoing drug development efforts aimed at intercepting and treating COVID-19, and help prepare for future coronavirus outbreaks. HIGHLIGHTSFocused CRISPR screens targeting host factors in the SARS-CoV-2 interactome were performed for SARS-CoV-2, HCoV-229E, HCoV-NL63, and HCoV-OC43 coronaviruses. Focused interactome CRISPR screens achieve higher resolution compared to genome-wide screens, leading to the identification of critical factors missed by the latter. Parallel CRISPR screens against multiple coronaviruses uncover host factors and pathways with pan-coronavirus and virus-specific functional roles. The number of host proteins that interact with a viral bait protein is not proportional to the number of functional interactors. Novel SARS-CoV-2 host factors are expressed in relevant cell types in the human airway.

Enhanced SARS-CoV-2 Neutralization by Secretory IgA in vitro (preprint)

SARS-CoV-2 primarily infects cells at mucosal surfaces. Serum neutralizing antibody responses are variable and generally low in individuals that suffer mild forms of the illness. Although potent IgG antibodies can neutralize the virus, less is known about secretory antibodies such as IgA that might impact the initial viral spread and transmissibility from the mucosa. Here we characterize the IgA response to SARS-CoV-2 in a cohort of 149 individuals. IgA responses in plasma generally correlate with IgG responses and clones of IgM, IgG and IgA producing B cells that are derived from common progenitors are evident. Plasma IgA monomers are 2-fold less potent than IgG equivalents. However, IgA dimers, the primary form in the nasopharynx, are on average 15 times more potent than IgA monomers. Thus, secretory IgA responses may be particularly valuable for protection against SARS-CoV-2 and for vaccine efficacy.

Escape from neutralizing antibodies by SARS-CoV-2 spike protein variants (preprint)

Neutralizing antibodies elicited by prior infection or vaccination are likely to be key for future protection of individuals and populations against SARS-CoV-2. Moreover, passively administered antibodies are among the most promising therapeutic and prophylactic anti-SARS-CoV-2 agents. However, the degree to which SARS-CoV-2 will adapt to evade neutralizing antibodies is unclear. Using a recombinant chimeric VSV/SARS-CoV-2 reporter virus, we show that functional SARS-CoV-2 S protein variants with mutations in the receptor binding domain (RBD) and N-terminal domain that confer resistance to monoclonal antibodies or convalescent plasma can be readily selected. Notably, SARS-CoV-2 S variants that resist commonly elicited neutralizing antibodies are now present at low frequencies in circulating SARS-CoV-2 populations. Finally, the emergence of antibody-resistant SARS-CoV-2 variants that might limit the therapeutic usefulness of monoclonal antibodies can be mitigated by the use of antibody combinations that target distinct neutralizing epitopes.

Measuring SARS-CoV-2 neutralizing antibody activity using pseudotyped and chimeric viruses (preprint)

The emergence of SARS-CoV-2 and the ensuing explosive epidemic of COVID19 disease has generated a need for assays to rapidly and conveniently measure the antiviral activity of SARS-CoV-2-specific antibodies. Here, we describe a collection of approaches based on SARS-CoV-2 spike-pseudotyped, single-cycle, replication-defective human immunodeficiency virus type-1 (HIV-1) and vesicular stomatitis virus (VSV), as well as a replication-competent VSV/SARS-CoV-2 chimeric virus. While each surrogate virus exhibited subtle differences in the sensitivity with which neutralizing activity was detected, the neutralizing activity of both convalescent plasma and human monoclonal antibodies measured using each virus correlated quantitatively with neutralizing activity measured using an authentic SARS-CoV-2 neutralization assay. The assays described herein are adaptable to high throughput and are useful tools in the evaluation of serologic immunity conferred by vaccination or prior SARS-CoV-2 infection, as well as the potency of convalescent plasma or human monoclonal antibodies.

Convergent Antibody Responses to SARS-CoV-2 Infection in Convalescent Individuals (preprint)

During the COVID-19 pandemic, SARS-CoV-2 infected millions of people and claimed hundreds of thousands of lives. Virus entry into cells depends on the receptor binding domain (RBD) of the SARS-CoV-2 spike protein (S). Although there is no vaccine, it is likely that antibodies will be essential for protection. However, little is known about the human antibody response to SARS-CoV-21-5. Here we report on 149 COVID-19 convalescent individuals. Plasmas collected an average of 39 days after the onset of symptoms had variable half-maximal neutralizing titers ranging from undetectable in 33% to below 1:1000 in 79%, while only 1% showed titers >1:5000. Antibody cloning revealed expanded clones of RBD-specific memory B cells expressing closely related antibodies in different individuals. Despite low plasma titers, antibodies to three distinct epitopes on RBD neutralized at half-maximal inhibitory concentrations (IC50s) as low as single digit ng/mL. Thus, most convalescent plasmas obtained from individuals who recover from COVID-19 do not contain high levels of neutralizing activity. Nevertheless, rare but recurring RBD-specific antibodies with potent antiviral activity were found in all individuals tested, suggesting that a vaccine designed to elicit such antibodies could be broadly effective.

LY6E impairs coronavirus fusion and confers immune control of viral disease (preprint)

Zoonotic coronaviruses (CoVs) are significant threats to global health, as exemplified by the recent emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)1. Host immune responses to CoV are complex and regulated in part through antiviral interferons. However, the interferon-stimulated gene products that inhibit CoV are not well characterized2. Here, we show that interferon-inducible lymphocyte antigen 6 complex, locus E (LY6E) potently restricts cellular infection by multiple CoVs, including SARS-CoV, SARS-CoV-2, and Middle East respiratory syndrome coronavirus (MERS-CoV). Mechanistic studies revealed that LY6E inhibits CoV entry into cells by interfering with spike protein-mediated membrane fusion. Importantly, mice lacking Ly6e in hematopoietic cells were highly susceptible to murine CoV infection. Exacerbated viral pathogenesis in Ly6e knockout mice was accompanied by loss of hepatic and splenic immune cells and reduction in global antiviral gene pathways. Accordingly, we found that Ly6e directly protects primary B cells and dendritic cells from murine CoV infection. Our results demonstrate that LY6E is a critical antiviral immune effector that controls CoV infection and pathogenesis. These findings advance our understanding of immune-mediated control of CoV in vitro and in vivo, knowledge that could help inform strategies to combat infection by emerging CoV.